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§ ■ ABSTRACT 

<N ' 

| We introduce a new statistic to measure more accurately the cosmic sound speed of 

Q ' clusters of galaxies at different redshifts. This statistic is evaluated by cross-correlating 

\Q • cosmic microwave background (CMB) fluctuations caused by the Sunyaev-Zel'dovich 

effect from observed clusters of galaxies with their redshifts. When clusters are dis- 
tributed in redshift bins of narrow width, one could measure the mean squared cluster 
peculiar velocity with an error cr C 2 <(300km/s) 2 . This can be done around z > 0.3 with 
| clusters of flux above 200 mJy which will be detected by PLANCK, coupled with high 

resolution microwave images to eliminate the cosmological part of the CMB fluctua- 
tions. The latter can be achieved with observations by the planned ALMA array or the 
NSF South Pole telescope and other surveys. By measuring the cosmic sound speed 
and the bulk flow in, e.g., 4 spheres of ~ lOO/i^Mpc at z = 0.3, we could have a direct 
measurement of the matter density 0.21 < Q m < 0.47 at 95% confidence level. 

6 , 

£j , Subject headings: Cosmic Microwave Background. Cosmology: theory. Cosmology: 

^ ! observations. 
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1. Introduction 

The Mach number A4 characterizes the coldness of the velocity field. It was originally intro- 
duced by Ostriker & Suto (1990) as a test of cosmological models. The cosmic Mach number is 
independent of the normalization of the mass density power spectrum and is insensitive to red- 
shift and bias in linear theory. It effectively measures the slope of the power spectrum up to the 
scales corresponding to the size of the region and can be used to constrain cosmological models 
for structure formation (Suto & Fujita 1990). Strauss et al. (1993) made a direct comparison of 
observations and models. They analyzed two samples of spiral and one of elliptical galaxies and 
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obtained M ~0.5-l, with results differing with galactic type. Nagamine et al. (2001) studied the 
dependence of the cosmic Mach number on overdensity, galaxy mass and age and concluded that 
either the Local Group was in a relatively low density region or the true mass density was Q, m ~0.2. 

Measurements of the velocity field using galaxies as distance indicators are affected by system- 
atic errors (Strauss & Willick, 1998) making difficult the determination of the Mach number on 
scales beyond 10 4 km/s. This can be avoided by analyzing measurements of the Cosmic Microwave 
Background (CMB) as proposed in this Letter. Hot gas in moving clusters produces redshift inde- 
pendent cosmic microwave background (CMB) anisotropies (Sunyaev & Zel'dovich, 1972 hereafter 
SZ). The thermal component (TSZ) is usually larger than the kinetic term (KSZ), but has a distinct 
spectral signature and can be removed using multifrequency observations, allowing the determina- 
tion of the cluster peculiar velocities (Holzapfel et al. 1997; Mauskopf et al. 2000; LaRoque et al. 
2002). At present, the systematic errors are still 3 to 5 times larger than the measured velocities. 
Recently, Benson et al. (2003), used SuZIE-II measurements of 6 clusters at z >0.2 to constrain the 
amplitude of their bulk motion, leading to the first such determination at intermediate redshifts. 

The Mach number is the ratio of the bulk flow component, which can be measured directly 
from the KSZ data alone (Haehnelt & Tegmark 1996, Kashlinsky k, Atrio-Barandela 2000, Aghanim 
et al. 2001), to the cosmic sound speed. We show that the cosmic sound speed can be reliably 
measured by cross-correlating the CMB temperature field at cluster locations with their redshifts. 
CMB temperature anisotropies and cluster redshifts encode information on peculiar velocities, but 
have different systematic errors. Combining the two quantities allows to measure the average sound 
speed of clusters. 



2. The Cosmic Mach Number and Cosmic Sound Speed. 

The cosmic Mach number is defined to be M(r,R) = |^ lk ; ^ . The bulk flow component is 
the mean squared velocity over a region of characteristic scale r and Csjiow is the sound speed in 
the reference frame of the mean flow. They are: V 2 ulk (r) = HQf 2 (Q m ,0,A)^ J °° P(k)W 2 (kr)dk, 
and C 2 = H 2 f 2 (n m ,n A )^ / °° P(k)W 2 (kR)dk, with C| flow = C 2 - V 2 n l Here f(n m ,n A ) ~ 
f2^ 6 + ^a(1 + ^m/2)/70 (Lahav et al. 1991), P(k) is peculiar mass power spectrum and W(x) is 
the window function of the survey. Cs depends mainly on the co-moving scales R of the collapsed 
moving objects. Clusters are regions of high density and their Cs may differ from the average. In 
the linear regime, Bardeen et al. (1986) showed that the rms peculiar velocity of peaks is smaller 
than that of field points, but the difference is less than 3% for R < 3/i _1 Mpc. Hereafter, we compute 
the rms peculiar velocity of clusters using linear theory with a Gaussian kernel and R = 1.5/t _1 Mpc. 
Borgani et al. (2000) verified that this window and scale give a good fit to the measured peculiar 
velocities of 18 clusters from the ENEAR catalog of peculiar velocities of elliptical galaxies. 

At present the error bars on peculiar velocities from KSZ (cr Vp ) are much larger than the esti- 
mated velocities vp, preventing reliable measurement of the Mach number using CMB data alone. 
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We show that measurements of the CMB temperature of clusters and of their redshift can be com- 
bined to determine M in some specific configurations. The temperature field at a cluster location 
includes (Carlstrom et al. 2002): the intrinsic CMB signal, foreground residuals, noise and TSZ 
(STxh = TksTx /m e c 2 ) and KSZ (STxin = rvp/c) components. In these expressions r the cluster 
optical depth to electron scattering, kpTx the intracluster electron temperature in eV, and vp the 
radial component of the peculiar velocity. All components have different amplitude depending on 
the beam of each experiment and, except the kinetic and intrinsic CMB temperature anisotropies, 
all have different frequency dependence. The TSZ component is the dominant contribution to the 
CMB signal at Z<3000 (<4') (Atrio-Barandela & Miicket 1999, Molnar k Birkinshaw 2000). Multi- 
frequency measurements at cluster locations remove the temperature anisotropy down to the KSZ 
signal plus some residual ST = Tqt(vp/c) + r; the latter is defined having zero mean and variance 
< r 2 >= a\ T . The uncertainty on the measured peculiar velocity of any individual cluster would 
be o~ Vp = c(cfat/T t). Presently, the measurements of the KSZ (clusters with r ~ 3 — 5 x 10~ 3 ), 
have uncertainty a Vp ~ lOOOkm/s (Benson et al. 2003, Carlstrom et al. 2002). 

The bulk flow velocity can be evaluated directly from the KSZ signal by adding the velocity 
of all clusters in the sample (Benson et al. 2003) with uncertainty: 



The rms peculiar velocity along the line of sight of clusters located on a thin shell in redshift 
space can be reliably measured by cross-correlating the temperature field at cluster locations with 
their measured redshifts. For this purpose we chose bins of width Az and for each cluster we define 
5z = z — z, where z is the average redshift of the clusters in the bin. Cluster redshifts can be 
determined by photometric/spectroscopic measurements of redshift of member galaxies or using 
SZ morphology (Diego et al. 2002). The latter would require only SZ observations but provide 
an accuracy of a z ~ 0.04(1 + z) for one object, while multiobject spectroscopic measurements of 
galaxies provide redshifts with uncertainties a z < 5 x 10~ 4 (Yee, Ellingson & Carlberg 1996). Thus 
redshifts of clusters in the prospective samples can be assumed to have no uncertainties. We also 
assume that clusters are homogeneously distributed within bins. Thus the cross-correlation gives: 



Hereafter, Vbuik an d Cg will be the radial components of the earlier quantities. 

Peculiar velocities of clusters and their distances are not independent variables. Clusters are 
chosen to be in shells of width Az in redshift space; this selection criteria introduces a correlation 
between d and vp: clusters that are most separated will have larger peculiar velocity if they are 
to be in the sample. If cluster velocities are typically vp ~ Cg, on average one would expect the 




(1) 



However, a Vp 3> Cg, so the error on the sound speed is a C 2 oc a 2 p and does not scale as iV 



{— ■ Sz) = (t) ± (t— ■ 5z) ± (r • 5z). 
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spread of clusters on real space to be: H Q Ad ~ cAz + 2Cs- Hence: 

(c5T-c8z) _ 2 cAz (cAz + 2C s ) 

T (r) ~ C ^ low ± '"'Jim {p) 7^ • (3) 
In this expression, a Vp is the average error on the measured velocity of members of the cluster 
sample and (v p ) = Vtmik if the region is smaller than the coherence length of the velocity field; 
otherwise (v p ) = Cs/yfNd- 

To test eq. (3) we analyzed a numerical simulation based on a P 3 M code with 256 3 particles 
(Faltenbacher et al. 2002) of a ACDM model with A = 0.7, a 8 = 0.87 and Hubble constant 
H Q = 70 km s _1 Mpc _1 . We used a "friends-of-friends" algorithm identifying clusters with (1.4 — 
23) x 10 14 /i~ 1 M o . We randomly selected 100 shells of width Az = 0.01 containing ~100 clusters 
each and found that, on average, less than 10% of the clusters were outside the distance interval 
cHg 1 . In 80% of the cases \ {tvp ■ cSz)\ < (t)(v p ) ( cAz J^f s ^ ; for the remaining 20%, the spread was 
a few percent larger. In all cases, the second noise term in eq. (3) was negligible. 

The error bar of Cs depends linearly on Az. To reduce it we compute Csjiow i n shells with a 
small spread in redshift. For this geometry, the dominant contribution is 



<^,„™ = I****./.; [motifs \Ju> 



Since bulk flows in shells are small, the sound speed in the shell and the matter reference frames 
are very similar: C| j low ~ C|. Compared with a direct estimation of C| from KSZ data alone, the 
cross correlation with redshift represents an improvement of more than one order of magnitude. 



3. Observational Prospects. 

For small errors (Barlow, 1989): {(t m / 1 M) 2 = (oy bulk /Vb u ik) 2 + (a C 2 /2C|) 2 . In order to obtain 
an accurate determination, bulk flow velocities have to be measured in regions much smaller than 
the sound speed. We propose to compute the Mach number by estimating the sound speed in a shell 
of width Az located at z s h e ii and to measure bulk flows in spheres of radius R centered at the same 
redshift. Observations of clusters to measure their SZ amplitude are already (almost) routinely 
carried out by several telescopes (e.g. Carlstrom et al. 2002) and exponential progress is expected 
with the construction of the ALMA array 4 and the NSF- funded South Pole microwave telescope 5 . 
For the purpose of this project, high resolution observations of large samples of clusters with well 
measured redshifts are required. Blind surveys searching for clusters in large areas of the sky, like 
the cluster catalog expected to be obtained from the full sky CMB observations of the upcoming 
PLANCK 6 mission or by the South Pole telescope, which will effectively observe ~ 10000 deg 2 , 



4 http://www. alma.nrao.edu/ 
5 http://astro. uchicago.edu/spt/ 

6 http://astro. estec.esa.nl/SA-general/Projects/Planck/ 
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will be most relevant. Clusters that produce a change in flux of about 200 mJy relative to the 
mean flux of the CMB will be detected by Planck. For the WMAP ACDM concordance model 
this limit translates into ~ 9000 clusters in 4/5 of the sky (Diego et al. 2003). Cluster redshifts 
could be coarsely determined using SZ data (Diego et al. 2002). Optical follow-up of those clusters 
without optical counterpart will be necessary to measure their redshifts. The current technologies 
have the capacity to measure ~2000 redshifts per night with an average uncertainty smaller than 
130 km/s 7 . If 10 galaxies per cluster are used to define the cluster redshift, this procedure would 
require a reasonable amount of observing time. 

Most clusters detected by PLANCK will be unresolved. Those needed for computing the 
Mach number will have to be re-observed at high resolution with ALMA, where its high fre- 
quency coverage and its small projected noise of ~ 2/xK rms after one hour of integration time 
for a beam of l', could be very useful to subtract foreground and CMB contribution to its min- 
imum. The flux limit can be related to the cluster mass (Kay et al. 2001): SV^SSGHz) = 



Planck channel, the one with the highest resolution and lowest noise. In this relation, M200 cor- 
responds to the mass within a region where the over-density is 200p c . Using the mass- luminosity 
relation (Reiprich & Bdhringer, 2002) L x ,boi = 1-15 x 1O 45 (M 2O o/1O 15 /i^ 1 M ) 1 - 8 /i~ 2 ergs s" 1 we 
can estimate the luminosity of clusters detected by Planck at any given redshift. 

Fig. 1 shows the relative errors on C| and A4 vs z. The number of clusters were modeled 
with the X-ray Luminosity Function (XLF) as an evolving Schechter function: (j)(L, z)=<p {l + 
z) L~ a exp(L/L*), with L*=L*^(\ + z) B , where A, B are two evolutionary parameters, and L*, a 
are the local XLF values. We took L*=9.3x 10 u h~ 2 ergs s" 1 and o =2.44x lO^/i^Mpc" 3 (Ebeling 
et al. 1997). In (a) and (b), solid lines corresponding to no redshift evolution (A, B) = (0,0), dashed 
to (-3,0) and dot-dashed to (-1,-2), obtained from the analysis of the ROSAT Deep Cluster Survey 
and the Extended Medium Sensitivity Survey (EMSS), respectively. At present, the evolution of the 
bright end of the XLF is not well determined and the degree of negative evolution deduced from the 
EMSS survey is most likely an overestimation, that does not extend to clusters with Lx < 10 44 /i -2 
ergs s" 1 (Rosati et al. 2002). For the 200 mJy limit (Fig. la) the RDCS evolution parameters 
predict about 14000 clusters out to z=l, and for the 400 mJy, about 6000. Those numbers were 
12000 and 4000 for the EMSS evolution parameters, respectively, the largest difference being in 
the number of clusters at z >0.3. The number density of the RDCS parameters was very similar 
to that of Kay et al. (2002). We took (7,^=1000 km/s, independent of z, even though foreground 
and CMB residuals and source confusion may depend on the cluster redshift. Any z dependence 
can be easily accounted for since both a C 2 and om vary linearly with a Vp , the average error on the 
velocity of individual clusters in the sample. 

Spheres of radius ~100/i _1 Mpc will contain very few clusters in blind all-sky surveys such as 



7 http: //www. sdss.org/, http://msowww.anu.edu.au/2dFGRS/ 




This detection threshold is quoted for the 353 GHz 
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Planck cluster catalog, and another type of survey is required. Romer et al. (2001) argue that an 
XMM serendipitous cluster survey of about ~800 deg 2 will detect more than 8000 clusters ranging 
from poor to very rich systems. Their effective flux limit would be ~ 1.5 x 10 _14 ergs s _1 cm -2 , and 
all clusters above kTx > 4KeV (or luminosities Lx > 10 44 /i -2 ergs s _1 ) will be detected. Adopting 
the t-Lx relation from Cooray (1999) and Kashlinsky & Atrio-Barandela (2000), rescaled to the 
bolometric band, the number cluster density (n c ) verifies that r(n c i?^) 1 / 2 ~1 is almost independent 
of luminosity, so that ~100 clusters per sphere of radius 100/i _1 Mpc with r*2 x 10~ 2 will be 
detected, independent of z. Also, they will have, on average, a similar error a Vp . In Fig. lc and Id 
we plot the error on the Mach number assuming that the sound speed was computed in shells of 
Planck detected clusters and the bulk velocity is computed using a top hat window on spheres of 
100 (thick) and 150 /i _1 Mpc (thin lines) radii. Solid, dashed and dotted-dashed lines represent the 
same XLF evolution models as in Figs. la,b. Fig. lc corresponds to Planck detecting all clusters 
with flux above 200 mJy and (d) to 400 mJy. We also considered cr„ p =1000 km s _1 , although 
now the bulk velocity and sound speed are computed over different samples. Folding any redshift 
dependence of a Vp into our results is straightforward, but would require detailed knowledge of the 
(unknown) cluster selection function for Planck and XMM. 

To see the capacity of our method to determine cosmological parameters, we assume that at 
redshift z ~0.3 we have identified all PLANCK clusters with fluxes above 200 mJy in a shell of 
Az=0.01, have their z measured with uncertainty smaller than 150 km/s and have them re-observed 
with ALMA at 1' resolution. For the expected number of 100 clusters, this would require ~100 
observing hours. The sample is used to compute Cs at that z. We also assume that the clusters 
identified by XMM in a sphere of radius 100/i _1 Mpc centered at z have data of similar quality. The 
Mach number would be the ratio of the bulk velocity in the sphere to the sound speed in the shell. 
Fig. 2 plots M. vs matter density for flat models with cosmological constant, S7a + ^m=l> and for 
three different values for the primordial spectral index n = 0.95, 1, 1.3 shown as dashed, solid and 
dot-dashed lines (all other cosmological parameters were chosen from the WMAP measurements, 
Bennett et al. 2003, Spergel et al. 2003). The matter power spectrum was taken as P(k) oc k n T 2 (k) 
with the analytical approximation of Sugiyama (1995) for the CDM transfer function, T(k), and 
normalized to as=0.8. The upper set of lines corresponds to bulk flows on spheres of 100/i~ 1 Mpc 
and the lower set to 150 /i _1 Mpc radius. For f2 m =0.3 the triangle shows the expected Mach number 
and its 1-a error bars. We assumed that the cluster XLF evolves with RDCS parameters. The 
error bar on M translates into a confidence interval for Q m . If the bulk flow is measured for 4 
independent spheres of 100/i~ 1 Mpc radius, we can estimate 0.25(0.21) < £l m < 0.37(0.47) and at 
the 68 (95) % confidence level. Our method is insensitive to a running spectral index, quintessence 
or tilt: variations on these cosmological parameters lead to differences in the Mach number much 
smaller than om- In this respect, our measurement of the matter density would be robust. 
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4. Conclusions. 



We presented a new method to measure the cosmic sound speed of clusters of galaxies. Assum- 
ing that all clusters with SZ flux larger than 200 mJy will be identified by the Planck mission, and 
that the XRLF evolution in the Rosat Deep Cluster Survey is representative of the overall cluster 
evolution, the relative error would be Oqi /(7| ~ 0.3. If the XMM Serendipitous Cluster Survey 
detects all clusters above kTx = 4KeV, by combining measurements of bulk flows in spheres of 
different sizes with the cosmic sound speed measured in shells at the same redshift, we can estimate 
the Cosmic Mach number with relative accuracy o^jM. > 0.2. This accuracy can be achieved 
with high resolution microwave images in order to reduce the intrinsic CMB anisotropy with re- 
spect to the Kinematic SZ signal, and is within the reach of the currently planned South Pole and 
ALMA telescopes. The Cosmic Mach number we can determine by our method is directly related 
to the matter density, and it can determine 0.21 < Q m < 0.47 at the 95% confidence level. This 
result, being independent of other mass density measurements, would provide an important self 
consistency check. 
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HA2002-0084. FAB also thanks the Junta de Castilla y Leon research project SA002/03 and 
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Fig. 1. — (a) Relative error of the sound speed for three different evolutions of the luminosity 
function: no evolution (solid line), dashed to (A,B) = (—3,0) and dot-dashed to (—1,-2) -see 
text-. The detection limit is 200 mJy. (b) Same as before, but the threshold limit is 400 mJy. (c) 
Relative error on the Mach number, for the same evolution histories as before. In every pair, the 
upper line corresponds to bulk flows on a sphere of 100/t _1 Mpc, and the lower line to a sphere of 
150/t _1 Mpc. The flux limit is 200 mJy. (d) The same as in (c) but the clusters detected have fluxes 
above 400 mJy. 
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Fig. 2. — Mach number estimated by measuring bulk flows on spheres of 100 (upper set) or 
150 /i _1 Mpc (lower set). Dashed, solid and dot-dashed lines correspond to a spectral index n = 
0.95, 1.0 and 1.03, in agreement with WMAP results. The triangles show the prediction for the 
WMAP concordance model. The la error was computed assuming that there were 4 independent 
measurements of the bulk flow at that redshift. They are slightly shifted for a better display. 
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